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ABSTRACT 

Knowledge of the chemical composition and absolute masses of Capella are key to understanding 
the evolutionary state of this benchmark binary system comprising two giant stars. Previous efforts, 
including our own 2009 study, have largely failed to reach an acceptable agreement between the 
observations and current stellar evolution models, preventing us from assessing the status of the 
primary. Here we report a revision of the physical properties of the components incorporating recently 
published high-precision radial velocity measurements, and a new detailed chemical analysis providing 
abundances for more than 20 elements in both stars. We obtain highly precise (~0.3%) masses of 
2.5687 ± 0.0074 M e and 2.4828 ± 0.0067M©, radii of 11.98 ± O.57R 0 and 8.83 ± 0.33 R Q , effective 
temperatures of 4970 ± 50 K and 5730 ± 60 K, and independently measured luminosities based on the 
orbital parallax (78.7 ± 4.2 L 0 and 72.7 ± 3.6 Lq). We find an excellent match to stellar evolution 
models at the measured composition of [Fe/H] = —0.04 ± 0.06. Three different sets of models place 
the primary star firmly at the end of the core helium-burning phase (clump), while the secondary is 
known to be evolving rapidly across the Hertzprung gap. The measured lithium abundance, the C/N 
ratio, and the 12 C/ 13 C isotopic carbon abundance ratio, which change rapidly in the giant phase, 
are broadly in agreement with expectations from models. Predictions from tidal theory for the spin 
rates, spin-orbit alignment, and other properties do not fare as well, requiring a 40-fold increase in 
the efficiency of the dissipation mechanisms in order to match the observations. 

Subject headings: binaries: general - binaries: spectroscopic — stars: abundances — stars: evolution 
— stars: fundamental parameters — stars: individual (Capella) 


1. INTRODUCTION 

As one of the brightest binary stars in the sky, 
Capella (a Aurigae, HD 34029, HR 1708, HIP 24608, 
G8 III+G0 III, P or b = 104 days, V = 0.07) has been 
studied for more than a century with a wide range of 
techniques and at all observable wavelengths^ A per¬ 
sistent source of frustration for several decades has been 
the difficulty in determining accurate absolute masses for 
the components, despite the wealth of astrometric and 
spectroscopic measurements available. The history of 
this problem has been related by sever al authors (e.g., 
iBatten et all I199H iBarlow et al.l I1993I1. and most re¬ 
cently in our earlier paper (iTorres et al.l I2009L hereafter 
T09). The challenge associated with the masses has hin¬ 
dered efforts to pin down the precise evolutionary state 
of the more massive primary star, which has widely been 
considered to be a core helium-burning object, based 
mostly on timescale arguments. Disappointingly, current 
stellar evolution models have so far largely failed to con¬ 
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5 Capella has a wide common proper motion companion that 

is itself a visual binary composed of M dwarfs. The system is 

therefore a hierarchical quadruple. Revised properties of the M 

dwarf pair are reported in the Appendix. 


firm that notion, as it has not been possible to achieve 
a satisfactory fit to the global properties of both stars 
simultaneously at a single age when assuming the bulk 
chemical composition the system is believed to have. The 
secondary, on the other hand, is clearly on its way across 
the Hertzprung gap. 

The uncertainty in the masses was thought to have 
been solved in the T09 analysis, which improved the for¬ 
mal precision by about a factor of three compared to 
previous estimates, and documented efforts to control 
systematic errors in the radial velocities that have likely 
plagued the determination of the velocity amplitude of 
the rapidly-rotating secondary star for decades, as de¬ 
scribed there. Despite this, it was still not possible to 
establish the state of the primary component unambigu¬ 
ously when enfor cing a sin g le age. _ 

In the interim, IWeber fe Strassmeiefi (1201 111 have pre¬ 
sented a new spectroscopic study of Capella based on 
much higher-quality observational material that leads to 
significantly larger masses for both stars than in our 2009 
study, by several times the stated uncertainties. In par¬ 
ticular, the spread in mass between the stars increased 
from 1% to about 3.5%, which is an enormous differ¬ 
ence for a pair of giants, and could drastically change 
the assessment of their relative state of evolution. In ad¬ 
dition, we have now made a new determination of the 
chemical composition of Capella that is appreciably dif¬ 
ferent from the abundance assumed in the earlier pa¬ 
per, and is a key ingredient for the comparison with 
stellar evolution models. These two developments mo¬ 
tivate us to take a fresh look at the system in order to 
investigate the impact of the new measurements. Fur- 
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thermore, iWeber fc Strassmeierl (120111) have presented 
evidence that the orbit of Capella may be very slightly 
eccentric, unexpectedly, whereas all previous studies in¬ 
cluding our own assumed it is circular. It is of inter¬ 
est, therefore, to revisit our 2009 study of tidal evolution 
in the system (orbit circularization and rotational syn¬ 
chronization) with more sophisticated models than used 
there, especially in light of the new masses. 

We have organized the paper as follows. In Sect. [2] 

we d escribe t he_new spe ctroscopic observations of 

IWeber fc Strassmeierl 1 2011 ) and comment on the issue 
of systematics in the radial velocities, in comparison with 
our previous results from 2009. A revised orbital fit for 
Capella is the presented in Sect. [31 using also astromet¬ 
ric and other measurements from T09. Sect. [4] reports 
a new detailed chemical analysis of both stars from dis¬ 
entangled spectra, along with a determination of the at¬ 
mospheric parameters. The revised physical properties 
of the stars are collected in Sect. 0 and are compared 
against three different sets of stellar evolution models in 
Sect. [6] Key chemical diagnostics available for Capella 
are also compared with model predictions in this sec¬ 
tion. Then in Sect. [7] we examine the evolution of orbital 
and stellar properties subjected to the influence of tidal 
mechanisms, as a test of that theory. Sect. [8] presents a 
discussion of the results and concluding remarks. Finally, 
the Appendix provides an update on the orbital proper¬ 
ties of the wide common proper motion companion of 
Capella. 


2. RADIAL VELOCITIES 

The numerous historical radial-velocity (RV) measure¬ 
ments of Capella have been discussed at length in our 
T09 study, which highlighted how challenging it has been 
to determine accurate values for the rapidly rotating 
secondary star (usini Ri 35 kms -1 ), whereas those of 
the sharp-lined primary (usini « 4 kms -1 ) have been 
quite consistent over the last century. T09 presented 
162 new RV determinations for both components based 
on spectra obtained at the Harvard-Smithsonian Cen¬ 
ter for Astrophysics (CfA) covering only a very narrow 
wavelength range (45 A). The RVs were measured using 
the t wo-dimensional cross-co rrelation algorithm TOD- 
COR dZucker fc Maxell 1 1 9941) . with synthetic templates 
appropriate for each star. Because of the limited wave¬ 
length coverage, those measurements are susceptible to 
systematic errors resulting mostly from lines shifting in 
and out of the spectral window as a function of orbital 
phase. Therefore, an effort was made to control those 
biases by performing numerical simulations to determine 
corrections to the velocities, which were at the level of 
the final uncertainties in the individual measurements for 
the secondary, and slightly larger for the primary. Final 
errors in the RVs as measured from the scatter in the 
orbital fit were about 0.44 kms -1 for the primary and 
0.89 kms -1 for the secondary. A sign that systematic 
errors remained at some level in the CfA velocities was 
evident in the residuals of the secondary star shown in 
Figure 2 of T09, in which a pattern can be seen as a 
function of orbital phase, with a peak semi-amplitude of 
about twice the typical error. Possible explanations for 
this, as discussed by T09, include the presence of spots 


on the active secondary star, or template mismatch^ An 
additional indication of possible biases was the fact that 
a small offset (0.267 ± 0.079 kms -1 ) was found between 
the primary and secondary velocities in the global orbital 
fit of T09 that could not be accounted for by differences 
in the gravitational redshift between the stars, and was 
ascribed to similar reasons as the secondary residual pat¬ 
tern. 

More recently IWeber fe Strassmeierl (l201ll ) have re¬ 
ported new RV measurements for Capella based on a 
very large set of more than 400 s pectra obtained with the 
STELLA echelle spectrograph (|Strassmeier et alJ 1201 0) 
on a 1.2 m robotic telescope in Tenerife, Spain. These 
spectra are of far superior quality than those of T09, both 
in terms of wavelength coverage (t wo orders o f magnitude 
larger) and signal-to-noise ratios. IWeber fc Strassmeierl 
(|2011H derived velocities using a similar two-dimensional 
cross-correlation approach as T09, and also performed 
numerical simulations to assess and correct for system¬ 
atic errors caused by the measuring technique. The ve¬ 
locity scatter from their orbital fit is 0.064 kms -1 for 
the primary and 0.297 kms -1 for the secondary, seven 
and three times smaller than in T09, respectively. The 
key difference in this data set compared to T09 is in 
the resulting velocity semi-amplitude of the secondary 
star (Kb = 26.840 ± 0.024 kms -1 ), which is more 
than 6cr larger than reported by T09 (Kb = 26.260 ± 
0.087 kms -1 ). This difference alone leads to absolute 
masses for Capella that are 4% larger than in T09 for 
the primary, and 2% for the secondary, a very significant 
change that exceeds the formal mass errors by a factor of 
many. The semi-amplitudes Ka of the primary star, on 
the other hand, are in virtually perfect agreement (see 
below). 

Despite t he muc h imp roved random errors of 
IWeber fe Strassmeied (120111) . the residuals of the sec¬ 
ondary velocities from their spectroscopic orbital model 
(see their Figure 2) still display a phase-dependent pat¬ 
tern reminiscent of the one in T09, also with a semi¬ 
amplitude of roughly twice the errors, but at a much 
lower level in absolute terms. Moreover, they note that 
a small offset is seen again between the primary and 
secondary velocities (0.059 kms -1 ) that cannot be ex¬ 
plained by the gravitational redshift effect. This sug¬ 
gests that systematic errors may still be lurking in these 
new measurements, for possibly some of the same rea¬ 
sons as before. Nevertheless, any remaining systemat¬ 
ics are likely to be significantly less important than in 
T09, as expected not only from the much higher qual¬ 
ity of the spectroscopic material, but suggested also by 
the significantly smaller magnitu des of _7) t he correc¬ 
tions for systematics applied by IWeber fe Strassmeierl 
dlmD, ii ) the formal uncertainties in the individual RVs, 
or equivalently, the scatter from the spectroscopic or¬ 
bital solution, in) the amplitude of the residual pattern 
for the secondary, and iv ) the unexplained residual pri¬ 
mary/secondary offset. In the next section we therefore 
incorporate these measurements in a new orbital analysis 

6 In particular, due to limitations in the available library of syn¬ 
thetic spectra they used, the macroturbulent velocity of the tem¬ 
plates (Cut = 1-5 kms -1 ) was not quite as large as appropriate for 
giant stars. This also resulted in an overestimation of the rotational 
velocities of the components, as discussed by T09. 
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of Capella. 

Of the 438 R V measureme nts reported in Table 1 
of IWeber fe Strassmeierl (j2011h . their final solution ex¬ 
cluded 14 for the primary and 7 for the secondary. We 
have done the same here as we found them to give un¬ 
usually large residuals, and we adopted also the measure¬ 
ment uncertainties as published. 

3. REVISED ORBITAL FIT 

The global orbital solution in our T09 study of Capella 
combined all usable astrometric observations in the liter¬ 
ature with the CfA RVs for both stars, as well as radial 
velocities for the primary star from many of the his¬ 
torical data sets. The latter were carefully examined 
to ensure that they imply Aa values consistent with 
those from the CfA RVs in separate spectroscopic so¬ 
lutions using the same fixed orbital peri od (see T0 9, Ta- 
ble 2). A similar solution of the IWeber fc StrassmeieH 
(l 201 lh velocities shows that the primary semi-amplitude, 
A'a = 29.959 ± 0.005 kms -1 , is essentially the same 
as that from the CfA RVs, A'a = 29.96 ± 0.04 kms -1 . 
Therefore, for our revised glo bal so lution below we have 
incorporated the IWeber fe StrassmeieH (120111) measure¬ 
ments for both stars, the CfA velocities for the primary 
(but not the secondary), and the primary velocities from 
the same historical data s e ts as used in T09. Although 
the IWeber & StrassmeieH ( 20 IT ) observations certainly 
dominate by weight, the older measurements are still use¬ 
ful because they extend the baseline more than a century, 
improving the orbital period. 

The extensive astrometry available for Capella includes 
measurements made by many authors with the technique 
of lon g-ba seline interferometry, beginning with the work 
of IMerrilll (ll922ih as well as speckle interferometry, direct 
imaging, and the intermediate astrometric measurements 
from Hipparcos. These have all been described and tab¬ 
ulated in our earlier study, and we refer the reader to 
that work for details. So far as we are aware, no further 
astrometric observations have been published for Capella 
except for those of lHubv et al.1 (120131b which we do not 
use here, however, because of concerns expressed by these 
authors about possible systematic errors affecting their 
measurements. 

Our global orbital fit follows closely that described by 
T09, and includes the following parameters: orbital pe¬ 
riod (Porb), relative angular semi-major axis (a"), incli¬ 
nation angle (i 0 rb), eccentricity (e), longitude of perias- 
tron of the more massive and cooler star (waQ posi¬ 
tion angle of the ascending node for the equinox J2000.0 
(Tj 2 ooo)> time of periastron passage (T per i), center-of- 
rnass velocity ( 7 ), the velocity semi-amplitudes for each 
star (A'a and An), the angular semimajor axis of the 
photocenter (ap hot ), corrections to the Hipparcos catalog 
values of the sky position of the barycenter (Act cos <5, 
AS) at the mean catalog reference epoch of 1991.25, and 
corrections to the proper motion components (A p, a cos 5, 
A p$). To account for differences in the zero points of the 
various RV data sets r e lative to the primary star mea¬ 
surements of IWeber & StrassmeieH ( 201T . we have also 
solved for 10 velocity offsets, one for each set. An ad¬ 
ditional parameter, f p , was included to correct the scale 

7 Note that this is the fainter star at optical wavelengths (see 
T09), which we refer to as star ‘A’. 


of the angular separation measurements from two of the 
a strometric data sets (s ee T09 for details). 

IWeber fc StrassmeieH (| 201 lH discussed several adjust¬ 
ments made to their secondary velocities to place them 
on the same system as their primary velocities of Capella. 
These adjustments were intended to correct for differ¬ 
ences in the gravitational redshift of the two components, 
and other shifts of unknown origin (see above). From 
their discussion it is not entirely clear to us whether these 
constant shifts have been applied to the velocities they 
reported, so our global solution includes one additional 
offset, Aab, to account for possible residual effects. With 
this, the total number of adjustable parameters in our fit 
is 27. 

The solution in our T09 study assumed a circular or¬ 
bit for Capella, as have all previous analyses of the bi¬ 
nary. We noted, however, that there were hints of a non¬ 
zero eccentr icity in the interferometric measurements of 
iHummel et all (I1994D . though not in the CfA RVs or in 
any of the oth er data sets. We ascrib ed this to the trans¬ 
formation that IHummel et ahl (Il994i ) made between their 
original interferometric visibilities ( V 2 ) and the nightly 
relative positions in polar coordinates that they pub¬ 
lished, given that their own solution using the original 
visibiliti es indicated a circular orbit. As pointed out by 
IWeber fc StrassmeieH ( 201 l b however, a spectroscopic 
fit using their RVs also indicates a statistically signifi¬ 
cant non-zero eccentricity, of very nearly the same mag¬ 
nitude as we had seen, and even with a consistent longi¬ 
tude of periastron. This suggests we may have been too 
quick to dismiss the possibility of a non-circular orbit in 
T09, as unexpected as this may be for a pair of giants 
(one being a clump star) i n a period of 104 days (see, 
e.g.. iMassarotti et al.l [20 08). We discuss this further in 
Sect. □ For our new global solution we have chosen to 
allow the orbit to be eccentric, on the assumption that 
the effect is real. 

The parameters of our fit are presented in Table |I] 
along with other properties derived directly from the or¬ 
bital elements. With the exception of A'b and the quanti¬ 
ties that depend on it, the other results are rather similar 
to those in T09. Th e eccentri city is essent ially the same 
as derived bv IWeber fc Strassmeierl (120111 b 

4. SPECTROSCOPIC ANALYSIS 

Until recently the only detailed c hem ical analysis avail¬ 
able for Capella was that bv iMcWilliam ( 199 0). indicat¬ 
ing a sub-solar composi tion of [Fe/H] = —0.37 ± 0.22 
on the scale of iCrevesse (fl984l ). equivalent to [Fe/H] = 
— 0.20 ± 0.22 on a more modern scale in which the solar 
iron abundance is A(Fe) = 7.50@ This determination is 
presumably based on the sharp lines of the primary star, 
but there is no indication that the presence in the spec¬ 
trum of the nearly equally bright secondary was properly 
accounted for, and in addition, the analysis adopted an 
incorrect primary temperatures A new metallicity deter¬ 
mination was reported by iFuhrmaiml (l201ll) that gives a 
rather higher abundance of [ Fe/H] = + 0.05 ± 0.08, ap¬ 
parently on the scale of lGrevesse et all (Il996i ) in which 
the solar iron abundance is also A(Fe) = 7.50. This 

8 We use the standard abundance notation in which A(X) = 
log[n(X)/n(H)] + 12, where n(X) and n(H) are the numbers of 
atoms per unit volume of element X and of hydrogen. 
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TABLE 1 

Revised orbital solution for Capella. 


Parameter 

Value 

Adjusted quantities 

Po rb (days). 

104.02128 

± 0.00016 

a" (mas). 

56.442 

± 0.023 

*orb (deg). 

137.156 

± 0.046 

e . 

0.00089 

± 0.00011 

wa (deg). 

342.6 

± 9.0 

Hj 2000 (deg). 

40.522 

± 0.039 

T peri (HJD—2.400,000). 

48147.6 

± 2.6 

7 (kms -1 ) a . 

+29.9387 

± 0.0032 

A'a (kins -1 ). 

25.9611 

± 0.0044 

K& (km s -1 ). 

26.860 

± 0.017 

“phot ( mas ). 

2.14 

± 0.70 

Aa* (mas). 

-0.53 

± 0.81 

AS (mas). 

-0.37 

± 0.57 

A/d a cos 8 (mas yr — ■*-). 

+0.33 

± 1.00 

A gs ( mas y r_1 ). 

-0.04 

± 0.60 

f P b . 

1.0400 

± 0.0032 

Aab for WS (km s 1 ) c . 

+0.050 

± 0.013 

Ai (WS-C01) (kms" 1 ) ... 

-0.19 

± 0.14 

A 2 (WS-N00) (kms- 1 ) ... 

+2.35 

± 0.45 

A 3 (WS-G08) (kms" 1 ) ... 

-0.67 

± 0.26 

A 4 (WS-S39) (kms- 1 ).... 

-1.52 

± 0.16 

A 5 (WS-S53) (kms- 1 ).... 

+0.62 

± 0.14 

A 6 (WS-B86) (kms" 1 ) ... 

-0.10 

± 0.13 

A 7 (WS-S90) (kms- 1 ).... 

-2.07 

± 0.60 

Ag (WS-B91) (kms- 1 ) ... 

-0.62 

± 0.13 

A g (WS-B93) (kms- 1 ) ... 

+0.79 

± 0.10 

Aio (WS-T09) (kms- 1 ) ... 

+0.289 

± 0.035 

Derived quantities 

Ma (M 0 ) . 

2.5687 

± 0.0074 

M b (M 0 ) . 

2.4828 

± 0.0067 

q = M b /Ma . 

0.96653 

± 0.00062 

a (10® km). 

111.11 

± 0.10 

“ (Rq) . 

159.72 

± 0.15 

a (au). 

0.74272 

± 0.00069 

7r orb (mas). 

75.994 

± 0.089 

Distance (pc). 

13.159 

± 0.015 

p, a cos S (mas yr —1 ). 

+75.85 

± 1.00 

PS (mas yr- 1 ). 

-427.17 

± 0.60 

(WCOn/ . 

1.204 

± 0.060 


Note. — References for the RV offsets Ai to Aio are: C01 
= ICampbelH (flMIh: N 00 = INewalll <fl900l); G08 = I Good dl908h : 

S39 — Struve! d 19391): S53 — Istruve & Kilbvl d 19531); B 86 = 
iBeavers^^JCittcr (198j|); S90 == Shcherbakov ct al. (1990); B91 
— I Bat ter^ltHdTl 19911) ; B93 = [Barlo^^^aLrJ^93n ; and T09 = 
ITorres et al.| d2009l) . The physical c onsta nts used in the analysis 
are those specified by Torres ct al. ( 2 0101) . 

a On the reference frame of the RVs of Weber &: Strassmcier 

(f2onh . 

k Scale factor for the angular separation measurements bv IMerrilll 
d!922f) and IKulaginl d 19701) . 

c Primary/secondary offset for the IWeber &; StrassmeieH (|2011l) 
velocities (WS), in the sense (primary minus secondary). 

^ Flux ratio between the secondary and primary in the Hipparcos 
passband, derived from the angular semimajor axis, the semima¬ 
jor axis of the photocentric orbit as measured by the satellite, 
and the velocity semi-amplitudes (see T09). 

study is based on spectral synthesis applied to eight Fe I 
lines and one Fe II line from a single composite spectrum 
with some degree of line blending, using an unspecified 
primary/secondary flux ratio. 

Below we describe our new determination of the 
chemical composition and atmospheric parameters of 
Capella based on the t echnique of spectral disentangling 
(|Simon fe Sturmlll994D , which bypasses the line blending 
problems inherent in previous analyses that used com- 



Fig. 1.— Measured flux ratios for the components of Capella 
(cool star relative to hot star, i.e., primary relative to secondary). 
Values from T09 are indicated with dots and error boxes, in which 
the horizontal length of eac h box indicates the w avelength coverage. 
Other measurements from I Hu by et al.l 120131 1 are shown by the 
lighter gray squares in the range 0.611—0.843 /im. Overplotted is 
the pr edicted flux ratio based on synthetic spectra bv IHusser et all 
1 20131 1 scaled according to the radius ratio given by T09. 
posite spectra. 

4.1. Disentangling 

For our spectroscopic analysis we made use of pub¬ 
lic archival spectra take n in 2003 and 2004 with the 
ELODIE spectrograph (iBaranne et al.1 I1996T) on the 
1.93 m telescope at the Observatoire de Haute-Provence, 
in France. The nominal resolving power of the instru¬ 
ment is R = 42, 000, and the 15 spectra used span the 
approximate wavelength range 4000-6800 A, with signal- 
to-noise ratios ranging from about 130 to 560 per pixel 
at 5550 A. We have disen tangled the se spe ctra using the 
FDBinary program of llliiic et al.l (l2004f) . in th e same 
way as described recently bv ITorres et all (|2014fl . FD¬ 
Binary implements spectral disenta ngling in the Fourier 
domain according to lHadraval (|1995l ). The signal-to-noise 
ratios of the resulting disentangled spectra are approxi¬ 
mately 510 for the cooler primary and 590 for the sec¬ 
ondary. Renormalization of the disent angle d spectra f or 
a proper abundance analy sis fsee lPavlovski fe Hensbergel 
120051 : iLehman et aLll2013fl requires knowledge of the rel¬ 
ative flux contribution of each star at each wavelength. 
Flux ratios for Capella have been measured throughout 
the UV, optical, and near-infrared range, as reported by 
T09. Figure |T| shows the predi cted flux ratio based on 
PHOENIX model spectra from iHusser et al.l (12013! ) for 
parameters near those of the components, along with the 
measurements tabul ated b y T09 as w ell as others from 
the recent study by iHubv etahl (|2013l ) over the range 
6112-8430 A. The agreement is very good. For our pur¬ 
poses we have used a smoothed version of this relation. 
Sample segments of the disentangled spectra of the two 
components of Capella are presented in Figure [2] 

4.2. Atmospheric parameters and abundance analysis 

The general methodology for determining the atmo¬ 
spheric parameters and abundances from the disen¬ 
tangled spectra follows the procedures described by 
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Fig. 2.— Portions of the disentangled ELODIE spectra of 
Capella. Secondary shifted vertically for clarity. 

Torres et al.1 (|2014f) . The uclsyn code (|Smallev et al.1 
200 1 [) was used to synthesize spectra, based on initial 
temperatures and surface gravities from T09 an d a b uilt- 
in grid of LTE model atmospheres bv ICastelli fc Kuruczl 
H2003H with a scaled-solar mixture. Line broadening was 
modeled adopting an initial composition matching the 
Sun, and microturbulent velocities of = 1.5 km s ^ 1 for 
both components (fine-tuned below). Least-squares fit¬ 
ting then yielded a macroturbulent velocity for the pri¬ 
mary of Ca = 6.6 ± 0.4kms _1 , and a projected rota¬ 
tional velocity of i>Asin* = 4.4 ± 0.5kms _1 . Given the 
rapid rotation of the secondary, macroturbulence has a 
negligible effect on the line profiles and we made no at¬ 
tempt to determine it from the ELODIE spectra. The 
rotational velocity of the companion was measured to be 
i>b sin* = 34.5 ± 6.7kms -1 . 

Elemental abundances were determined using spec¬ 
tral lines suitable for giant stars over the wavelength 
range 4600-6750 A. Line_lists an d at omic data were taken 
from the wo rk oflReddv et ah (l2012fl iBocek Toncu et ahl 
uni), and iLvubimkov et all (120151 ) Equivalent widths 
measured within UCLSYN for the very numerous iron lines 
were used to set the microturbulent velocities from the 
condition of a null correlation between the abundance 
and the reduced equivalent widths. We derived val¬ 
ues of £t,A = 1-48 ± 0.08kms _1 for the primary and 
£t,B = 1.55 ± 0.11 kms -1 for the secondary. We also 
made an estimate of the effective temperatures from 
the usual condition of excitation equilibrium, iterating 
with the measurement of £ t , with the following results: 
T eff , A = 4980 ± 80 K and T eStB = 5750 ± 110 K. There 
is very good agreement between these values and others 
reported by T09; we discuss them further below. The 
surface gravities in our analysis were held fixed at the 
estimates reported by T09, which are very close to our 
final values described in the next section. 

Detailed abundances were determined for 22 species in 
both stars, and oxygen in the primary only. They are 
listed in Table [2j which incl udes also the va l ues re lative 
to the Sun on the scale of lAsplund et al.1 (12009D . No 
adjustments have been applied for NLTE effects. The 
uncertainties account for possible errors in T e g as re¬ 
ported above, and also include a contribution from a rep¬ 



Atomic number 

Fig. 3.— Photospheric abundance pattern measured for the 
Capell a com p onents, c ompared to the standard solar composition 
of lAsplund et ahl 112009 ) (gray shading). Abundances from different 
ions of the same element have been averaged. 

resentative error of 0.1 kms^ 1 in £ t - The uncertainties in 
log g have a negligible impact. The choice of the mixture 
adopted in the model atmospheres, particularly the CNO 
composition, also has a minimal effec t on our abundance 
determinations (see also lMorel et all 120141 ). We find no 
dependence of the Fe abundance on wavelength, which 
is an indication that our adopted wavelength-dependent 
flux ratios (Figure 1) are accurate and do not introduce 
significant systematic errors in the abundances. A fur¬ 
ther indication of the robustness of our determinations is 
the fact that the abundances are quite similar for the two 
components (except for species affected by evolution; see 
below), as expected for a binary system. 

The weighted average iron abundance of the two stars 
from Fe I is [Fe/H] = —0.04 ± 0.06, or very nearly so¬ 
lar, in contrast with the sub-solar composition adopted 
by T09, and in better agreement with the estimate by 
iFuhrmaun d201lh . We find no significant enhancement 
of the a elements in Capella: [a/Fe] = —0.02 ± 0.04. 
A graphical representation of the abundance pattern for 
the two components is seen in Figure |3l compared to the 
solar composition. 

The lithium abundance has long been known to be 
very different for the t wo components of this binary 
dWallersteinl 11964 I1966D as a result of chemical evolu¬ 
tion in th e primary. We used spe ctral lines and atomic 
data from ILvubimkov et al.l (|2012l ) in the vicinity of the 
Li I A6708 doublet to make new estimates for each 
star, and obtained values of A(Li) = 1.08 ±0.11 and 
A(Li) = 3.28 ± 0.13 for the primary and secondary, 
respectively. These are consistent with previous mea¬ 
surements. The equivalent widths we determined are 
21.4 ± 0.7 mA and 297.4 ± 8.0 mA. 

Additional chemical diagnostics for Capella have been 
reported by T09, and include the 12 C / 13 C carb on isotope 
ratio for the primary star (27±4: lTomkin et al.l 197(1 ) and 
the C/N abundance ratios for both components (0.57 ± 
0.06 for the primary and 3.30 ± 0.16 for the secondary). 
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TABLE 2 

Abundances from our disentangled ELODIE spectra of Capella. 


A 

X 

Primary 


Secondary 


l°ge 0 

Abundance 

[X/H] 

N 

Abundance 

[X/H] 

N 

3 

Li I 

1.08 i 0.11 

-0.07 i 0.15 

3 

3.28 + 0.13 

+2.30 + 0.16 

3 

1.05 + 0.10 

6 

C I 

8.25 i 0.14 

-0.18 i 0.15 

4 

8.28 + 0.11 

-0.15 + 0.12 

5 

8.43 ± 0.05 

8 

O I 

8.55 i 0.11 

-0.14 i 0.12 

1 




8.69 + 0.05 

11 

Nal 

6.13 i 0.09 

-0.11 i 0.10 

4 

6.33 + 0.07 

+0.09 + 0.08 

4 

6.24 ± 0.04 

12 

Mg I 

7.60 i 0.09 

i0.00i0.10 

2 

7.42 + 0.10 

-0.18 + 0.11 

5 

7.60 + 0.04 

14 

Si I 

7.69 i 0.04 

+0.18 i 0.05 

10 

7.59 + 0.07 

+0.08 + 0.08 

7 

7.51 + 0.03 

20 

Cal 

6.27 i 0.11 

-0.07 i 0.12 

8 

6.38 + 0.11 

+0.04 + 0.08 

7 

6.34 + 0.04 

21 

Sc I 

3.13 i 0.15 

-0.02 i 0.16 

5 

3.16 + 0.10 

+0.01 + 0.11 

6 

3.15 + 0.04 

21 

Sc II 

3.12 i 0.07 

-0.03 i 0.08 

8 

3.10 + 0.08 

-0.05 + 0.09 

11 

3.15 + 0.04 

22 

Ti I 

4.96 i 0.12 

+0.01 + 0.13 

15 

5.02 + 0.09 

+0.07 + 0.10 

11 

4.95 + 0.05 

22 

Ti II 

4.93 i 0.05 

-0.02 i 0.07 

3 

4.91 + 0.09 

-0.04 + 0.10 

6 

4.95 + 0.05 

23 

VI 

4.07 i 0.12 

+0.14 + 0.14 

14 

4.10 + 0.07 

+0.17 + 0.11 

13 

3.93 + 0.08 

24 

CrI 

5.64 i 0.09 

+0.00 + 0.10 

9 

5.67 + 0.07 

+0.03 + 0.08 

11 

5.64 ± 0.04 

24 

Cr II 

5.61 i 0.09 

-0.03 + 0.10 

7 

5.57 + 0.07 

-0.07 + 0.08 

6 

5.64 + 0.04 

25 

Mn I 

5.32 i 0.09 

-0.11 + 0.05 

8 

5.31 + 0.08 

-0.12 + 0.09 

5 

5.43 + 0.05 

26 

Fe I 

7.47 i 0.06 

-0.03 + 0.07 

42 

7.44 ± 0.08 

-0.06 + 0.09 

41 

7.50 + 0.04 

26 

Fe II 

7.39 i 0.07 

-0.11 + 0.08 

8 

7.38 + 0.06 

-0.12 + 0.07 

11 

7.50 + 0.04 

27 

Co I 

4.87 i 0.08 

-0.12 + 0.10 

8 

5.03 + 0.07 

+0.04 + 0.10 

5 

4.99 + 0.07 

28 

Ni I 

6.20 i 0.04 

-0.02 + 0.06 

16 

6.21 + 0.07 

-0.01 + 0.08 

17 

6.22 + 0.04 

39 

Y II 

2.11 i 0.09 

-0.10 + 0.10 

4 

2.23 + 0.05 

+0.02 + 0.07 

5 

2.21 + 0.05 

40 

Zr I 

2.54 i 0.07 

-0.04 + 0.08 

5 

2.25 + 0.12 

-0.33 + 0.13 

3 

2.58 + 0.04 

57 

La II 

1.11 i 0.08 

+0.01 + 0.09 

5 

1.23 + 0.05 

+0.13 + 0.06 

5 

1.10 + 0.04 

60 

Nd II 

1.49 i 0.06 

+0.07 + 0.07 

8 

1.52 + 0.05 

+0.10 + 0.06 

7 

1.42 + 0.04 


Note. — Columns list the atomic number, the element and ionization degree, the logarithm of the number 
abundance on the usual scale in which A(H) = 12, the logarithmic abundance relative to the Sun, and the number 
of spectral lines measured. The last column gives the reference photospheric solar values from [Asplund et al.| 
f2009l) . 


We adopt those as published. 

5. PHYSICAL PROPERTIES 


TABLE 3 

Revised physical parameters of Capella. 


The revi sed masses of C ape lla that incorporate the 
new RVs of lWeber fc Strassmeierl ( 2011 ) have formal un¬ 
certainties of 0.3%, and differ only slightly from theirs 
through a combination of a marginally larger Kq value in 
our analysis and a smaller orbital inclination angle than 
they used. The new masses are 4.2% and 1.6% larger 
than those listed by T09, which is a significant difference 
due almost entirely to the change in the velocity ampli¬ 
tude of the secondary. The radii, based on the individual 
angular diameters from T09 and the revised orbital par¬ 
allax, are approximately 0.9% larger than before, and 
have precisions of about 5% and 4% for the primary and 
secondary, respectively. These are limited by the angu¬ 
lar diameters, as the uncertainty in the orbital parallax 
is only 0.11%. 

T09 reported three independent estimates of the effec¬ 
tive temperatures for the two components. One is from 
a comparison of their spectroscopic observations with 
synthetic spectra with solar metallicity, giving T e s t a = 
4900 ± 100 K and T e g ,b = 5710 ± 100 K. Another came 
from the use of the measured angular diameters of the 
stars along with their apparent magnitudes, the parallax, 
and bolometric corrections. The updated parallax in the 
present work does not alter those estimates significantly; 
they are T^a = 4970 ± 160 K and T^b = 5690 ± 130 K. 
A third determination by T09 was based on the mea¬ 
sured color indices for the stars, and th e use o f the 
color/temperature calibrations of Ramirez fe Melendez! 
(120051) . which depend on metallicity. A sub-solar com¬ 
position [m/H] = —0.37 ± 0.07 had been assumed by 
T09, whereas we now derive a value much closer to so¬ 
lar. Using our determination of [Fe/H] = —0.04 ± 0.06 


Parameter Primary Secondary 


Mass (M 0 ). 2.5687 ± 0.0074 2.4828 ± 0.0067 

q = Mb/M a . 0.96653 ± 0.00062 

a (10 6 km). 111.11 ± 0.10 

a (au). 0.74272 ± 0.00069 

7r or b (mas). 75.994 ± 0.089 

Distance (pc). 13.159 ± 0.015 

Radius (i? 0 ). 11.98 ± 0.57 8.83 ± 0.33 

logg (cgs). 2.691 ± 0.041 2.941 ± 0.032 

T eff (K). 4970 ± 50 5730 ± 60 

Luminosity (L 0 ) a ... 78.7 ± 4.2 72.7 ± 3.6 

BC V (mag). -0.304 ± 0.055 -0.089 ± 0.051 

M v (mag). 0.296 ± 0.016 0.167 ± 0.015 

v sin i (km s —1 ) b . 4.1 ± 0.4 35.0 ± 0.5 

P rot (days) c . 104 ± 3 8.5 ± 0.2 

Age (Myr) d . 590-650 

[Fe/H]. -0.04 ± 0.06 

A(Li). 1.08T0.11 3.28 i 0.13 

12 C/ 13 C d . 27 i 4 

C/N c . 0.57 i 0.06 3.30 i 0.16 


a Computed from V. 7T 0 rb ■ and BCv from iFloweil tf 1996 4 ■ adopting 
M® ol = 4.732 (see T09 and lTorrejl201Cli . 

d Average of 5 measurements from the literature for the primary and 
10 for the secondary that account for macroturbulence, including our 
own (see text). 

c Measured values adopted from T09. 

^ Age range from the MESA and Granada models (see text). 
e Measurement bv lTomkin et all J1976D . 


(Sect. 14. 21) . the revised photometric estimates become 
4940 K and 5680 K, which are 30 K and 70 K higher than 
before. Furtherm ore, a car eful ex amina tion of the zero 
point of the iRamfrez fe Melende j (12005D calibrations by 
iCasaerande et al.l ( 20101 ) suggests that the scale of those 
relations is too cool by about 85 K compared to the best 
available absolute scale, at least in the temperature range 
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of the Capella components. We have therefore applied 
this offset, obtaining corrected photometric estimates of 
T eff , A = 5025 ± 110 K and T eff)B = 5765 ± 120 K. The 
uncertainties include a contribution of 100 K added in 
quadrature to the photometric and calibration errors, to 
be conservative. Finally, a fourth temperature determi¬ 
nation was reported in the previous section from the dis¬ 
entangled ELODIE spectra, giving T e ff ; a = 4980 ± 80 K 
and T e ff,B = 5750 ± 110 K. The weighted average of the 
four values for each component is T e ff a = 4970 ± 50 K 
and T e ff b = 5730 ± 60 K, in which the uncertainties ac¬ 
count not only for the individual weights but also for 
the scatter of the measurements, and are believed to be 
realistic. These averages are 50 K hotter than in T09. 

The available determinations of v sin i for both com¬ 
ponents were summarized in our previous work (T09, 
Table 14). Our present measurements from the ELODIE 
spectra are consistent with those of others, as well as 
with the measurements reported by iFuhrmannl (120 1 ill , 
which are UAsini = 3.5 ± 0.8kms _1 and ussini = 
35.4 ± 3.2 km s . The weighted averages of all inde¬ 
pendent determinations (5 for the primary, 10 for the 
secondary) that have taken account of macroturbulence 
broadening, especially for the primary component, are 
va sin i = 4.1±0.4kms _1 and vb sin * = 35.0±0.5kms _1 , 
which we adopt for the remainder of the paper. 

The masses, radii, temperatures, and other derived 
properties are summarized in Table [3] Note that the 
bolometric luminosities are independent of the temper¬ 
atures and radii, and are based on the apparent mag- 
nitudes, the o rbital parallax, and bolometric corrections 
from iFlowerl (11996D . as in T09. If we instead compute 
them from T e g and R, the results are consistent, but 
have larger formal uncertainties: La = 78.7±8.1 Lq and 
Lb = 75.4 ± 6.4 L 0 . 

6. COMPARISON WITH STELLAR EVOLUTION MODELS 

Up until now the ability of stellar evolution models to 
match all of the global properties of both components 
of Capella simultaneously at a single age has not been 
entirely satisfactory, likely at least in part because there 
are so many observational constraints available. This 
has made it difficult to establish the evolutionary sta¬ 
tus of the primary star unambiguously, although it has 
widely been thought to be a core helium-burning (clump) 
star, based on timescale arguments (see T09). The sig¬ 
nificantly different (and more precise) masses obtained 
above, and evidence that the chemical composition is 
rather different from that previously assumed, motivate 
us to revisit the comparison with stellar evolution models 
here. 

An initial test was ma de using the PARSEC isochrones 
of iBressan et al.l ( 2012 ) FI These models adopt the solar 
distribution of hea v y elem ents from the compilation by 
iCrevesse fc Sauval (il998l). with adjustments to some el- 
ements following iCaffau et al.l (120111) such that the solar 
photospheric metallicity is Zq = 0.01524. On this scale 
the measured abundance of Capella (Sect. 14.21) corre¬ 
sponds approximately to Z = 0.0133. The helium abun¬ 
dance Y follows an adopted enrichment law with a slope 
AY/ AZ = 1.78, which results in a value for Capella of 
Y = 0.272. Convection is treated in the standard mix- 

9 http://stev.oapd.inaf.it/cmd . 


ing length theory approximation. The calibration to the 
Sun leads to a mixing length parameter of «mlt = 1.74, 
which is held fixed in these models. Convective core over¬ 
shooting is also included, with an efficiency parameter of 
A c = 0.5 for Capella, representing the mean free path 
of convective bubbles across the border of the convective 
region, expressed in units of the pressure scale height H p . 
This is roughly equivalent to a ov = 0.25 pressure scale 
heights above the convective boundary in the more com¬ 
monly used formulation of this phenomenon (see below). 
Mass loss from stellar winds is not considered in these 
models for the mass range of interest for Capella, nor is 
rotation. 

With the chemical composition and convective param¬ 
eters fixed as described above, we searched for the com¬ 
mon age giving the best fit to the masses, radii, tem¬ 
peratures, and luminosities of both stars using the y 2 
statistic 


X 



'ATeffl 2 

. CT Te ff . 



<?L 


’ Ai?l 2 

. °R _ 


as the figure of merit, where the sum is over both stars 
and the A quantities represent the difference between 
the predicted and measured properties. A reasonably 
good fit was obtained for an age of 622 Myr, matching 
the properties of the stars within about 1.4 times their 
uncertainties, with the largest discrepancy being in the 
primary temperature. We note, however, that the mass 
ratio of Capella is known much more precisely than the 
individual masses, with cr q ~ 0.06% compared to mass 
errors om of 0.29% and 0.27%. The masses in the above 
fit were allowed to vary independently, and as a result 
the best-fit mass ratio differs from the measured value 
by about 3er. We therefore repeated the fit constraining 
q to be near its measured value by using a corresponding 
penalty term [A q/a q ] 2 in y 2 instead of the secondary 
mass term. We obtained a solution of similar quality (all 
properties reproduced within 1.4cr) and nearly the same 
age of 625 Myr. This fit is illustrated in Figure [H and it 
places the primary star in the core helium-burning phase 
(clump). 

Aside from the global structural properties of the stars 
considered above, several chemical diagnostics including 
the lithium abundance, the C/N ratios, and the isotopic 
carbon abundance ratio 12 C/ 13 C are available for Capella 
that are not normally tabulated with published stellar 
evolution models, but that are nevertheless interesting 
to compare with theoretical predictions. To that end, we 
have perfor med an additional test against a second set of 
models by iClaretl (I2004D . occasionally referred to below 
as the Granada models. 

These models adopt also t h e sola r abundance distri¬ 
bution by IGrevesse fe Sauvall (11998D . with some adjust¬ 
ments such that the solar metallicity becomes Zq = 
0.0189. The abundance of Capella then corresponds ap¬ 
proximately to Z = 0.0172, which we held fixed. The 
enrichment law AY/AZ = 2.0 typically used in these 
models results in a helium abundance for Capella of Y = 
0.274, similar to that used above. The solar-calibrated 
value of the mixing length parameter is Omlt = 1-68, and 
convective core overshooting a ov is parametrized such 
that the mean free path above the convective boundary 
is d ov = a ov H p . Rotation was initially not included in 
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Fig. 4.— Comparison of the observed properties of Capella (M, 
R , T e ff, L, q) with a PARSEC isochrone from the model series by 
IBressan et al.l (120121) in the L— T e g and R-T e g- planes. Metallic- 
ity and age are as indicated in the lower panel. The insets show 
enlargements around the primary star, with a short line segment 
connecting the measured location in each plane with the best-fit po¬ 
sition along the isochrone. Corresponding insets for the secondary 
are not shown as the match to its properties is better. 

our tests. 

A grid of Granada evolutionary tracks was computed 
for the measured masses and a range of convective pa¬ 
rameters for each component of Capella, as in principle 
there is no reason to expect stars in such different evo¬ 
lutionary states to have the same convective properties. 
We varied ctmlt between values of 1.0 and 2.2 in steps 
of 0.1, and a ov over the range 0.15-0.40, with a step of 
0.05. Mass loss was inclu ded in these calculations fol¬ 
lowing the prescription by iReimersI (|1975t l . Preliminary 
tests indicated a minimal loss of mass for both stars up 
to the present age, but we nevertheless incremented the 
initial values slightly by 0.005 M 0 and 0.002 M 0 , respec¬ 
tively, so as to reproduce the measured masses exactly 
at the best-fit age. An excellent fit to the radii, tem¬ 
peratures, and luminosities of both stars was found for 
mixing length parameters of 1.80 for the primary and 
1.50 for the secondary, and convective core overshooting 
parameters of 0.35 and 0.30, respectively, with estimated 
uncertainties in each of these of about 0.05. Deviations 
in R, T e g, and L from the measured values are all smaller 
than 0.4cr. The best-fit age we obtained, about 649 Myr, 
was constrained to be the same for the two stars. We 
illustrate this solution in Figures [5] and [6] for the L-T e g 
and R-T e ff diagrams, respectively. 

We point out that the age in this best-fit solution is 
driven entirely by the properties of the secondary, specif¬ 
ically, by its effective temperature. This is because that 



Fig. 5.— Comparison of the observed propertie s of Capella (d ots 
and error boxes) with the Granada models by I Claret! (120041) in 
the L-T e ff diagram. The evolutionary tracks shown are for the 
measured masses (incremented by 0.005 Mq for the primary and 
0.002 Mq for the secondary, to account for mass loss; see text) and 
the measured metallicity (Z = 0.0172 for these models). Reference 
ages are marked along the secondary track, and the insets show 
enlargements around the position of each star, with a short line 
segment connecting the observations to the predicted positions. 

star is in such a rapid phase of evolution that the temper- 
ature is predicted to change drastically (by many times 
the observational uncertainty) in just 1 Myr. Figures [5] 
and [6] show, for example, that between the ages of 648 
and 649 Myr the temperature of the secondary decreases 
by 1130K, cooling by a further 660 K over the next mil¬ 
lion years. The primary, on the other hand, stays at es¬ 
sentially the same temperature over this time, and only 
changes its radius and luminosity, but at a much slower 
pace (see top insets in the figures). Consequently, it does 
not constrain the age nearly as much. Because of the 
rapid evolution of the secondary, the formal uncertainty 
in the age that comes from its temperature error is neg¬ 
ligible. A more meaningful uncertainty may be obtained 
by varying its mass within allowed limits, which results 
in an age range of approximately ±5 Myr. This does not 
include possible systematic errors having to do with the 
physics in the models. As found above from the PAR¬ 
SEC models, the primary star is seen to be located in 
the clump, on the hot side of the giant loop (end of core 
helium-burning phase), where the radius and luminosity 
are increasing with time. 

A test with Granada models that include rotation for 
both components gave a very similar fit, with an age of 
655 Myr that is only marginally older than before. 

The evolution of the surface chemistry in giants is the 
result of mixing and is directly related to the depth of the 
convection zone, which changes drastically as the stars 
approach the so-called first dredge-up, during their ini¬ 
tial ascent of the giant branch. Other changes can occur 
later. The first dredge-up event is illustrated in Figure [3 
which shows predictions from the Granada models for the 
abundance of lithium, the 12 C/ 13 C carbon isotope ratio, 
and the C/N ratio as a function of time. Also shown 
for reference are the changes in the location of the bot- 
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Effective Temperature (K) 


Fig. 6.— Similar to Figure[5]for the R-T e g plane. 

tom of the convection zone for each star (lower panel). 
The measurements of these key chemical diagnostics are 
represented with dots at the best-fit age of 649 Myr. 
Generally there is good agreement between theory and 
observation, except for the C/N ratios that deviate the 
mostH We note, however, that the C/N measurements 
rely on emission fluxes from spectral lines in the lower 
transition layers between the stellar chromosphere and 
the corona (specifically, C IV A1550.8 and N V A1238.8; 
see T09), so they may not strictly represent the abun¬ 
dances in the pho tosp h ere (despite some evidenc e that 
they do; see, e.g.. iBdhm-Vitense &; Mena-Werthl 119921) . 
The convection zone in the secondary star is seen to have 
just begun deepening, and should reach maximum depth 
approximately 7 Myr from now, according to these mod¬ 
els. 

A final test was performed against stellar evolu¬ 
tion tracks computed using the Modules for Experi¬ 
ments in Stellar Astrophysics (MESA, revision 7385; 
see IPaxton et all [MTT1 [ 2 ?Tl 3 ) FI Th ese models use the 
scaled solar abundances of lAsplund et all (12009H . accord¬ 
ing to which Z 0 = 0 .0134. The measured composition 
of Capella corresponds to a metal mass fraction of about 
Z = 0.012, and the helium abundance follows from an 
adopted /S.Y//S.Z = 1.67, and is Y = 0.2 70. Mas s loss 
is again computed according to the iReimer 1 (IT975h pre¬ 
scription, in this case with an efficiency parameter (mul¬ 
tiplicative scale factor) of ? 7 r = 0.2. For this paper 
we used t he ‘g rey_and_kap’ surface boundary condition 
(IPaxton et al.ll2013h . with o pacities and equation of st ate 
as discussed extensively bv IPaxton et al . (1201 4. 12013 ). 

Overshoot mixing across convective boundaries is 
treated slightly differently than in the models considered 
previously . MES A use s the exponential decay formalism 
of lHerwig et al.l (I1997T) . in which the product of the free 
parameter / ov and the local pressure scale height provide 
the scale length over which the diffusion coefficient de- 

10 The predicted difference between the ratios, however, is in 
better agreement with the measured difference (to within 1.7cr). 

11 http://mesa.sourceforge.net/ . 



ondary. The best-fit age of about 649 Myr is marked by the vertical 
dotted lines. The lower panel shows the depth of the bottom of 
the convective zone as a function of time, in units of the stellar 
radius. The dots represent the measurements, with error bars that 
are barely visible on this scale. Measured and predicted values 
are listed near each point, with the deviations listed in parenthesis 
next to the predictions, in units of the observational errors. 

cays from its value in the convective region. Although 
MESA allows for the free parameter to take on differ¬ 
ent values depending on the nuclear burning present in 
each convective zone, we have opted to use the same / ov 
value for H- and He-burning regions. For reference, with 
the choices listed above, we obtained a so lar-c alib rated 
mixing length parameter q:mlt = 1-84, and lHerwig et al] 
(jl997tl suggest a value of f ov ~ 0.02 for overshoot mixing 
above the H-burning core. 

A grid of MESA evolutionary tracks with the speci¬ 
fied composition was computed for each star over broad 
ranges in q:mlt (1-70 to 2.00, in steps of 0.05) and f ov 
(0.00-0.04, in steps of 0.01). The models were evolved 
from the fully-convective pre-main sequence to the end 
of core He burning. Due to the inclusion of mass loss, we 
increased the initial masses of the stars, as we did before 
for the Granada models, by 0.0048 M 0 and 0.0027 M 0 in 
this case such that the tracks reach the measured masses 
at their respective present locations in the H-R diagram. 
An excellent fit to the properties of both components was 
achieved at a common age of 588.5 Myr, with all residu- 
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Effective Temperature (K) 

Fig. 8. — Similar to Figure [5] for the MESA models. The mea¬ 
sured metallicity corresponds to Z = 0.0120 in these models. 



als being smaller than 1.2er. The optimal convective pa¬ 
rameters were found to be ctmlt = 1-85 for the primary 
and ctmlt = 1-80 for the secondary, with / ov = 0.02 for 
both stars. The matches in the L-T e g- and R-T e g dia¬ 
grams are illustrated in Figures [8] and [9] Once again the 
models place the primary at the end of the core helium¬ 
burning phase. The age is somewhat younger than ob¬ 
tained from the Granada models (a ^10% difference), 
while the age found earlier from the PARSEC models 
is intermediate between these two. These age differences 
correlate strongly and have to do mostly with the Z value 
used in the calculations for each model. The differences 
in Z at a fixed (measured) [Fe/H] value are in turn a con¬ 
sequence of the different heavy-elenrent mixtures adopted 
for the Sun in each case. 

The predictions from the MESA models regarding the 
evolution of the surface chemistry of the stars are similar 


to those from the Granada models. In particular, the 
predicted 12 C/ 13 C ratio for the primary at the age of 
588.5 Myr is 20.7, and the expected C/N ratios for the 
primary and secondary are 0.77 and 3.98, respectively, 
both being somewhat higher than measured. 

7. TIDAL EVOLUTION 

The considerable amount of information available for 
Capella offers a valuable opportunity to test our under¬ 
standing of tidal theory in binary stars, in ways those 
models have not often been challenged before. Such 
calculations are capable of making detailed predictions 
about the evolution of the size and shape of the orbit, as 
well as the rotational properties of the individual compo¬ 
nents including their spin rate, f2 rot = 27r/P rot , and the 
spin-orbit angle (f> (angle between the spin axis and the 
total angular momentum vector of the orbit, sometimes 
referred to as ‘obliquity’). In addition to the known or¬ 
bital elements of Capella, estimates are available also of 
the rotation periods of both stars, P ro t (see Tabled and 
T09), and of their projected rotational velocities, usin*. 
Our earlier study of the binary examined its tidal evo¬ 
lution considering the turbu lent di ssipat ion and radia¬ 
tive damping mechanisms bv [Zahnl (|1992t and references 
therein), as well as the hydrodynamical mechanism of 
iTassoul fe Tassoul (|1997l and references therein). These 
theoretical formulations involve a number of assump tions 
and simplifications discussed by 1 Zahnl ( 1977 ) and iHutl 
( 198i (h In particular, the equations are linearized around 
the equilibrium state, and are strictly valid only for rela¬ 
tively small eccentricities and near-synchronous rotation, 
as well as relatively small obliquities. 

In this work we have chosen to use the more gen¬ 
eral equations of tidal evolution by IHutl (119811) . which 
are valid for arbitrary eccentricities and rotation rates, 
although they are still restricted to relatively low mu¬ 
tual inclination angles <j). We used a fourth-order 
Runge-Kutta algorithm to integrate the six coupled 
differential equations describing the time-dependent 
changes in the orbital semimajor axis ( da/dt ), eccen¬ 
tricity ( de/dt ), angular rotational velocities of both 
stars (dflrot.A/dt, dfI ro t,B/dt), and their spin-orbit an¬ 
gles ( d(j)\/dt , d<j)B/dt). In what follows we normalize 
the rotation rates to the mean orbital rate 12 0 rbi for 
convenience. The relevant stellar properties that also 
evolve with time, such as the radius, were taken at each 
integration step d irect ly from the best-fit evolutionary 
tracks of iGlaretl ( 2004 ) discussed in the previous sec¬ 
tion. The turbulent dissipation timescale for the stel¬ 
lar phases with convective enve lop es ( later stages for 
Capella) was adopted from I Zahnl (|1977l) . whereas the 
times cale for earlier p h ases w ith radiative envelopes fol¬ 
lows iClaret fc Cunhal ( 199 7). The initial conditions, 
which are unknown, were set by trial and error to match 
the measured values of the orbital period, the eccentric¬ 
ity, and the spin rates at the current age as closely as 
possible. 

The outcome of these calculations is illustrated in the 
top four panels of Figure flOl Setting the initial values 
to Po = 220 days, eg = 0.70, (fl r ot,A/^°rb)o = 320, and 
(f^rot,B/f^orb)o = 260 leads to evolved properties that are 
very close to those observed at the present age. In partic¬ 
ular, the observed super-synchronous rotation rate of the 
secondary (fI ro t,B/^orb = -Porb/-Prot,B ~ 12) is well re- 
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Fig. 10. — Predicted tidal evolution for Capella according to 
inns cemd, compared with the observations (filled circles; error 
bars are smaller than the point size), (a) Eccentricity evolution. 
The vertical dotted line in this and the other panels marks the 
current age of the binary (649 Myr) according to the models by 
E2 aretl (120041) described in the text, (b) Evolution of the orbital 
period, (c) and (d) Evolution of the spin rate of the primary and 
secondary, normalized for convenience in terms of the orbital an¬ 
gular velocity. The dot-dashed line represents the evolution of the 
pseudo-syn chronous rotation rate, computed following Eq.(42) of 
IHuU d 19811) . (e) Evolution of the fractional gyration radius of each 
star, for reference (solid line for primary, dashed for secondary), 
(f) Evolution of the spin-orbit inclination angle for each star (lines 
as in previous panel), (g) and (h) Predicted projected rotational 
velocities of the two stars. 

produced. Qualitatively the largest difference is perhaps 
in the orbital eccentricity, which theory predicts should 
strictly have fallen to zero some 35 Myr ago, driven al¬ 
most exclusively by the evolution of the primary star. 
Quantitatively, however, the difference in e between the¬ 
ory and observation is small, as the measured value (if 
real) is only e = 0.00089 ± 0.00011. While the agreement 
reached in these four observed properties is not entirely 
unexpected because we have allowed for four free parame¬ 
ters (the initial values), we note that the good fit was only 
possible by increasing the nominal efficiency of the tidal 
dissipation by more than an order of magnitude. With¬ 
out this increase, we find that the predicted rotational 
velocities of the stars near the zero-age main-sequence 
are unreasonably low (vsini ~ 20kms -1 ) for early A- 
type stars, such as the Capella components would have 
been. In order to yield more reasonable projected ro¬ 
tational velocities in excess of lOOkms -1 we had to in¬ 
crease the efficiency of the tidal mechanisms by a factor 
of ^40. A similar shor tcoming in the efficienc y of the¬ 
ory was found earlier bv lClaret fc Cunhal (119971) . in their 
analysis of tidal synchronization and circularization of a 
sample of detached eclipsing binaries. 

Figure [TUI displays the evolution of the spin-orbit an¬ 


gles for the two stars near the present age, where the 
integrations have been performed with arbitrary initial 
values of (f >o = 0.4 radians (about 23°) for both stars, 
as we have no direct handle on those angles. Tests with 
other (non-zero) values show that the behavior is quali¬ 
tatively always the same: the primary’s spin axis aligns 
itself with the orbital axis (i.e., cf> reaches zero) well before 
the current age, whereas the secondary remains formally 
misaligned. The alignment of the primary happens at 
very nearly the same time that its rotation becomes syn¬ 
chronous and that the orbit circularizes. The spin-orbit 
angle of the secondary quickly shrinks to zero shortly af¬ 
ter the present age (4-5 Myr later), at which time it also 
synchronizes its rotation with the mean orbital motion. 
In both cases the changes are the result of significant 
structural adjustments in the stars, such as a sharp in¬ 
crease in the radius of gyration, r gyr , related to the mo¬ 
ment of inertia through I = M(r gyl R) 2 (see Figure [TUbl. 

While we cannot measure the present-day values of </>, 
it is possible to gain indirect knowledge about these an¬ 
gles using our spectroscopic estimates of v sin i along with 
the measurements of P ro t and the radii of both stars. 
These quantities are trivially related by 

• ■ 2itR ■ ■ fi\ 

vsmi = — —sinirot, (1) 

rot 

where the inclination angle on the left-hand side (inacces¬ 
sible to direct observation) is strictly that of the stellar 
spin axis relative to the line of sight (* ro t), and also ap¬ 
pears on the right. We may thus solve for the sini rot 
term on the right-hand side. Figure fill displays the dis¬ 
tributions of sin i ro t derived from the propagation of all 
observational errors in a Monte Carlo fashion. For ref¬ 
erence we also show the sine of the orbital inclination 
angle relative to the line of sight (i or b), which is directly 
measurable from the astrometric observations (Table [3]). 
The close agreement between sin ? TOt and sin i or b for both 
stars is highly suggestive that the spin axes may actually 
be parallel to the orbital axis in space, and this in turn 
would imply </> = 0. We cannot rule out a difference 
in quadrants such that the spin axes are tilted with re¬ 
spect to the axis of the orbit while still maintaining the 
same sine value (e.g., * ro t = 180° — « or b)> but such a co¬ 
incidence for both stars seems rather unlikelyQ This 
indirect empirical evidence that the obliquity may cur¬ 
rently be zero for both stars (which would hardly happen 
by chance) appears to point to a discrepancy with the 
prediction from tidal theory for the secondary compo¬ 
nent (Figure [TOf), whereas for the primary there is good 
agreement. 

It is possible that part of the reason for this differ¬ 
ence lies in the small-angle approximation implicit in the 
tidal differential equations of IHu 3 (i 1 98 li l for the angle 
(j>. However, we note also that the secondary of Capella 
happens to be in such a rapid state of evolution that the¬ 
oretical predictions for the spin-orbit angle are extremely 
sensitive to other details, including those related to struc¬ 
tural changes in the stars happening at this stage. Those 

12 Barring this type of situation, conversion of the distributions 
in Figurc lTTl to inclination angles yields i rotj a = 135.3^tg'g deg and 
*rot,B = 138.2^ 2 ’fi deg, compared to the orbital value of i OT \, = 
137.156 ±0.046 deg. 
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Fig. 11.— The histograms represent the empirical distributions 
of sinirot derived from the measurement of usini, P T ot, and the 
radius of each star as given in Table [3] The dashed lines mark 
the value of sin 2 or b determined from our orbital fit. The good 
agreement strongly suggests that the spin axes of both stars may 
in fact be parallel to the orbital axis. 

changes may actually have a larger impact than the dissi- 
pation processes themselves (though the latter must also 
matter). In particular, changes in both the spin rate and 
the spin-orbit angle are strongly driven in part by the 
sudden change in the moment of inertia (or equivalently, 
the gyration radius) illustrated in Figure ITOb . An ad¬ 
ditional complication is the large ad-hoc increase in the 
efficiency of the tidal mechanism that was required to 
match other observations, as mentioned earlier. Thus, a 
definitive assessment of the accuracy of tidal theory re¬ 
lated to its other approximations is difficult in this case. 

An alternate way of comparing observations with the 
predictions from tidal and evolution models combined is 
by examining the evolution of usini, which is a spec¬ 
troscopically measured quantity. It is a function of the 
stellar radius, the spin rate, and the inclination angle of 
the rotation axis to the line of sight, i ro t (Eq.pQ). The 
latter angle is not directly predicted by theory, but is 
related to other angles by 

cos 4> = cos i OI b cos i ro t + sin i OI b sin * rot cos A , (2) 

in which the obliquity <f> can be predicted, i OI b is known, 
and A is the angle between the sky projections of the 
spin axis and the orbital axis (‘projected obliquity’). The 
angle A changes with time and is directly measurable in 
eclipsing systems by observing the Rossiter-McLauglilin 
effect. As Capella does not eclipse we have no knowledge 
of this angle in this case, except for the weak condition 
that it represents a lower limit to the thre e-dimensional 
angle ^ (see, e.g.. iFabrvckv fe Winnll2009ll . Eq.[2] may 
be solved for the quantity sini ro t that we need in order 
to compute v sin i with eq. |T] , resulting in the quadratic 
equation 

A sin 2 i rot + B sin i Iot + C = 0 

where 

A = cos 2 i or b + sin 2 i OI b cos 2 A 
B = — 2 cos (j> sin * or b cos A 
C = COS 2 (j) — COS 2 iorb ■ 

For the primary star the prediction that (f> = 0 near the 


current epoch means that A is also zero. Using this value 
in the equations above to compute the expected evolu¬ 
tion of the projected rotational velocity leads to the trend 
shown in Figure ITTlb . where we refer to usini more prop¬ 
erly as vx sini rot ,A- Values at times when <j> ^ 0 will be 
somewhat less accurate because the projected obliquity 
A may also be different from zero. The predicted value 
of 4.0 kms -1 at the current age is in excellent agreement 
with the measurement (4.1 ± 0.4 kms -1 ). The corre¬ 
sponding evolution of WBsin* rot ,B is seen in Figure fTOlr. 
In this case the predicted value at the current age does 
not match the measurement, and the discrepancy is in 
part a reflection of the evolution of <f> that was discussed 
above, but also has to do with the very rapid changes in 
the structure of the star ( R , r gyr ) at the present time. 
Tests in which we changed A within reason yielded very 
similar results, and cannot explain the difference. 

Finally, it is of interest to verify that the components 
of Capella have always been detached, as any significant 
mass transfer earlier in their lives (e.g., through Roche 
lobe overflow in the primary) would invalidate our com¬ 
parison with stellar evolution models, which are designed 
for normal (unperturbed) stars. The fraction of its Roche 
lobe filled by the primary component depends on its size 
and also on the size of the orbit, which changes due 
to tidal forces (Figure [TUbl. The same Granada evo¬ 
lutionary track used above indicates the star attained a 
maximum radius at the tip of the giant branch of about 
38 Rq at an age of 617 Myr (MESA models predict a 
similar maximum size of 36 Rq), and the filling factor 
(R/Rt toche) at the time was approximately 0.63. This in¬ 
dicates that mass transfer through Roche-lobe overflow 
has not taken place in Capella. 

8. DISCUSSION AND CONCLUDING REMARKS 

For the first time current stellar evolution models are 
shown here to provide a satisfactory fit to the observed 
global properties of both components of Capella at a sin¬ 
gle age, for a chemical composition equal to that mea¬ 
sured. The comparison confirms the long-held but largely 
unsubstantiated belief that the primary is a clump star, 
and more precisely, it suggests Capella A is near the end 
of its core He-burning phase, as indicated consistently 
by three different sets of models. The principal fac¬ 
tors that have allowed the better match are an improve¬ 
ment in the accuracy of the absolute masses, made pos¬ 
sible by the high-qualit y radi al-velocity measurements 
of lWeber fc Strassmeierl ( 2011 ). and the first robust de¬ 
termination of the metallicity of Capella derived here. 
Our detailed chemical analysis of the disentangled spec¬ 
tra of the components yields essentially the same near¬ 
solar composition for the two stars, which is rather dif¬ 
ferent from the sub-solar abundance the binary was pre¬ 
viously assumed to have in our earlier T09 study. All 
measured properties (masses, radii, temperatures, and 
independently derived luminosities) are now simultane¬ 
ously in agreement with the models to within 0.4-1.4<r, 
depending on the model. This result strengthens our 
confidence in evolutionary calculations for evolved stars. 

Chemical indicators of evolution that differ greatly be¬ 
tween the components, such as the 12 C/ 13 C ratio, the 
lithium abundance, and the C/N ratio, are also broadly 
in agreement with theoretical predictions, though with 
somewhat larger differences that may be due to either 
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shortcomings in the mixing prescriptions in the models, 
or observational errors in these delicate measurements. 
Rarely has it been possible to perform this type of test 
involving key chemical diagnostics for a binary with two 
evolved components having well-determined absolute di¬ 
mensions (M, R, etc.), such as Capella. 

On the other hand, the performance of tidal theory 
as de scribed by the tidal evolution equations of iHutl 
lfl98l is not as good . As was noticed previously by 
[Claret fc Cunhal (1 99it ). we again find that the efficiency 
of the tidal mechanisms involved seems to be too low 
by a factor of ~40, if predictions are to be consistent 
with the large rotational velocities Capella A and B 
presumably had earlier in their lives as A-type stars 
(usin i > 100 kms -1 ), and at the same time the much 
lower velocities they now have as evolved stars. Appli¬ 
cation of an ad-hoc adjustment of this magnitude brings 
agreement, although theory still forecasts a significant 
misalignment between the spin axis of the secondary and 
the orbital axis at the present age, whereas empirical ev¬ 
idence based on the direct measurement of the star’s ra¬ 
dius, rotation period, and ysin* seems to favor an obliq¬ 
uity consistent with zero. Model predictions about the 
synchronous rotation of the primary and its spin-orbit 
alignment do agree with observational clues, as does the 
supersynchronous rotation of the secondary. A small ad¬ 
ditional disagreement in that the orbit seems to be very 
slightly eccentric even though it is expected to be per¬ 
fectly circular may or may not be significant, as the dif¬ 
ference is small. 

Finally, another important consequence of our revised 
chemical composition for Capella concerns the relation¬ 
ship with the coronal abundances, which have been mea¬ 
sured by many authors as summarized by T09. Coronal 
abundances in the Sun are known to depend on the first 
ionization potential (FIP), in such a way that low-FIP 
elements (less than 10 eV) are overabun dant compared 
to those with higher FIP (see, e.g., iFeldman fc Widina 


|2002T> . T09 showed this to be the case for Capella as 
well. Other stars display the oppo site effect (see, e.g., 
iBrinkman et al.l |2001k lLamind 120041) . What is less clear 
from evidence in other stars is whether the low-FIP el¬ 
ements are enhanced relative to the photospheric abun¬ 
dance, or whether it is the high-FIP elements that are de¬ 
pleted compared to the photosphere. T09 adopted a sub¬ 
solar photospheric composition for Capella of [m/H[ = 
—0.34 ± 0.07, based on the work of iMcWilliaml (119901 ). 
and since this agreed with the measured coronal abun¬ 
dances of the high-FIP elements (see their Figure 18), 
they concluded it was the low-FIP elements that were 
enhanced. Interestingly, our revised photospheric abun¬ 
dance of [Fe/H] = —0.04 ± 0.06 leads to precisely the 
opposite conclusion: now it is the low-FIP elements that 
agree with the photosphere, and therefore the high -FIP 
eleme nts are depleted. A recent study by iPeretz et all 
(IMl) of the coronae of stars with super-solar pho- 
tospheric abundances found the same effect shown by 
Capella in a Cen A and B, which have somewhat similar 
temperatures, although these are dwarfs so it is unclear 
how significant this may be. In any case, Capella now 
represents a robust point of reference for such studies in 
evolved and active stars. 

We thank Brian Mason (U.S. Naval Observatory) for 
providing the historical positional measurements of the 
Capella HL system, and the anonymous referee for help¬ 
ful comments. The research has made use of the SIM- 
BAD database, operated at CDS, Strasbourg, France, 
of NASA’s Astrophysics Data System Abstract Service, 
of the Washington Double Star Catalog maintained at 
the U.S. Naval Observatory, and of data products from 
the Two Micron All Sky Survey (2MASS), which is a 
joint project of the University of Massachusetts and the 
Infrared Processing and Analysis Center/California In¬ 
stitute of Technology, funded by NASA and the NSF. 


APPENDIX 

CAPELLA AS A MEMBER OF A MULTIPLE SYSTEM 

Over the last two centuries at least half a dozen visual companions to Capella hav e been recorded at angular 
separations ranging from 47" to 12'. Only the widest one, discovered by IFuruhieiml (I1914D . has been shown to 
be physically associated, and has received the double-star designation Capella H (also ADS 3841H, GJ 195 A, and 
LTT 11622, among others). The projected linear separation is ~9500 au (723" in position angle 141°). This 10th- 
magnitude star is of spectral type M1-M2.5, depending on th e source, and has the same relativel y large proper motion 
as Capella itself (I Furuhieiml 11914 iLepine fc Bqngior noll2007j ) as well as the same parallax (e.g., iDaniel fc Schlesinged 
1192ft I Adams et all 119261 : iBagnuolo fc Hartkopj 119891 ) . The earliest reported radial velocity measurements appear to 
be those of lAbtl (|1970l) at Mount Wilson, of which the first two from consecutive nights in 1921 gave a mean of 
about +30 kms -1 , while a third from 1927 is +48.5 kms -1 , with a large probable error. iStauffer fc Hartmannl ( 1986 ) 
reported a mean radial velocity of +30 ± 1 kms -1 based on an unspecified number of CfA spectra using the same 
setup as T09, and lTJpgren fc Carusol (|1988l) listed three measurements also fro m CfA taken in 1986, giving similar 
values. The first two of these ar e likely to be base d on the same spectra us ed by Stauffe r’fc Hartmannl (119861 ). Other 
measurements were reported bv lReid et alJ (119951 ) (+31.9 ± 15 kms - x ) and lGizis et all (120021 ) (+32.7 ± 1.5 kms 1 ). 

In order to supplement these values we report eight additional observations of Capella H with the same instrumen¬ 
tation and procedures described by T09. Radial velocities were derived with the IRA iFl task xcsao, using a synthetic 
template selected t o match the properties of th e sta r with T P ff = 3750 K and no rotational broadening. The three 
spectra obtained by IStauffer fc Hartmannl (I1986D and lUpgren fc Carusol (I1988D were reanalyzed in the same way, and 
our velocities supersede those in the original papers. We list all these m easurements in Table [4] on the absolute helio¬ 
centric frame defined by the IAU Radial-Velocity Standard Stars (see iStefanik et all! 1999 ). The 11 velocities show no 


13 IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which is operated by the Association of Universities 


for Research in Astronomy, Inc., under contract with the National 
Science Foundation. 
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TABLE 4 

New radial velocity measurements of Capella H. 


HJD 

(2,400,000+) 

Julian Year 

RV 

(kms -1 ) 

CRV 

(kms -1 ) 

46475.6230 .... 

1986.121 

+31.33 

0.82 

46494.5453 .... 

1986.173 

+30.97 

0.82 

46788.8338 .... 

1986.978 

+31.84 

0.70 

52986.8030 .... 

2003.947 

+32.15 

0.63 

53030.7455 .... 

2004.068 

+32.19 

0.65 

53055.6708 .... 

2004.136 

+32.16 

0.55 

53078.6029 .... 

2004.199 

+31.71 

0.51 

53106.5530 .... 

2004.275 

+31.43 

0.67 

53361.7995 .... 

2004.974 

+31.93 

0.58 

53419.6511 .... 

2005.133 

+31.09 

0.58 

53689.8944 .... 

2005.872 

+31.08 

0.56 



Fig. 12.— Measurements of the relative position of the Capella HL binary, along with the 388-yr orbit b ylfleintd (1 1 9751 ) ( dashed), and ou r 
revise d orbit (solid). Filled symbols represent measures from the Washington Double Star Catalog. IStearnsI 111933 1. and IHelminiak et all 
f 20091 ). and open symbols are for the photographic measurements by Heintz. The line of nodes in the revised orbit is indicated with a 
dotted line, and the motion on the plane of the sky is direct (arrow). The primary (Capella H) is represented by the plus sign. Short line 
segments connect the measurements with their predicted location on the revised orbit. 


variation within the measurement errors, and have a mean of +31.63 ± 0.14 kms' 1 . This is 1.69 ± 0.14 kms -1 larger 
than the center-of-mass velocity of Capella reported in Table [1] and is most likely due to motion in the ^9500 au 
orbit 0 

Capella H was found bv iSteariii (1193611 to have a close (~2") companion approximately 3.5 magnitudes fainteiFI. 
designated Capella L (also ADS 3841L and GJ 195 B). Significant changes in position relative to Capella H since the 
discovery date indicate orbital motion. Thus, the Capella system is a hierarchical quadruple. Historical measurements 
of the Capella HL binary from the Washington Double Star Catalog were kindly provided by Brian Mason, and are 
shown in Figure fl2l They were obtained mostly with visual micrometers or photographica lly. We have suppl emented 
these observations with additional photographic measurements reported b v IStea r nsI (ll939f) and iHeiiit^ (|1975f) . as well 
as mor e recen t and accurate adaptive optics measurements bv IHelminiak et al.1 ( 20091 ) . A tentative 388-yr orbit by 
iHeintd (Il975f) based on his photographic data and earlier observations is also shown, but does not seem to represent 
the bulk of the observations very well. 

While the astrometric orbit is still largely undetermined because of the short arc covered by the measurements, a 
helpful constraint may be obtained by making use of rough estimates of the masses of the stars d erive d from existing 
near-infrared photometry (2MASS) and the empirical mass-luminosity relation of lDelfosse et al.l (I2000H . With inferred 


14 Part of the difference (~0.5 kms -1 ) is due to the larger gravi¬ 
tational redshift of the M dwarf compared to the pair of giants, and 
there could also be a contribution from the difference in convective 
blucshifts. Additiona lly, there may be a sligh t offset between the 
zero-point of the IWeber fc Strassmeier 1 20111 1 velocities and the 
IAU system. 


15 This visual brightness estimate bv IKuiperl (1 1 93+ may be over¬ 
estimated, however, since the 2MASS JHK a magnitudes of the 
individual components differ by only a few tenths of a magnitude. 
Alternatively, the companion may itself be multiple. 
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values of approximately 0.57 M 0 and 0.53 M 0 for Capella H and L, and the use of the orbital parallax from Table [T] 
Kepler’s Third Law then leads to a constraint on the ratio a 3 /P 2 . A tentative new orbital solution derived in this 
way, which is also shown in Figure IT2l provides a significantly better fit with standard visual elements P 0 rb ~ 300 yr, 
a ~ 3"5, e ~ 0.75, i ~ 52°, w l ~ 88°, flj 2 ooo ~ 288°, and periastron passage near the year 2220. 
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